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ABSTRACT 
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Context. Data on the cosmic microwave background (CMB) radiation by the Wilkinson Microwave Anisotropy Probe 
(WMAP) had a profound impact on the understanding of a variety of physical processes in the early phases of the 
Universe and on the estimation of the cosmological parameters. 

Aims. Here, the 7-year WMAP data are used to trace the disk and the halo of the nearby giant spiral galaxy M31. 
Methods. We analyzed the temperature excess in three WMAP bands (W, V, and Q) by dividing the region of the 
sky around M31 into several concentric circular areas. We studied the robustness of the detected temperature excess 
by considering 500 random control fields in the real WMAP maps and simulating 500 sky maps from the best-fitted 
cosmological parameters. By comparing the obtained temperature contrast profiles with the real ones towards the M31 
galaxy, we find that the temperature asymmetry in the M31 disk is fairly robust, while the effect in the halo is weaker. 
Results. An asymmetry in the mean microwave temperature in the M31 disk along the direction of the M31 rotation 
is observed with a temperature contrast up to ~ 130 /iK/pixel. We also find a temperature asymmetry in the M31 
halo, which is much weaker than for the disk, up to a galactocentric distance of about 10° (~ 120 kpc) with a peak 
temperature contrast of about 40 /iK/pixel. 

Conclusions. Although the confidence level of the signal is not high, if estimated purely statistically, which could be 
expected due to the weakness of the effect, the geometrical structure of the temperature asymmetry points towards a 
definite effect modulated by the rotation of the M31 halo. This result might open a new way to probe these relatively 
less studied galactic objects using high-accuracy CMB measurements, such as those with the Planck satellite or planned 
balloon-based experiments, which could prove or disprove our conclusions. 
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, 1. Introduction Naturally, the importance of understanding the galactic 
Q\ • halos is closely related to the nature and distribution of 
0> Galax y rotation, in particular for the Andromeda galaxy the dark matter) which is re levant for the formation and 
i (M31) has been well studied especially in the optical, IR, dynamics of galaxies. In this respect, the methodology 
. and radio bands, and it gives important information on the adopte d in the present paper of using WMAP data to probe 
L; ■ mass distribution not only in galactic disks but also in their both the disk and the halo of M31j even if with the limi _ 
, halos UBiflBsd UaaaP. On the other hand, since they are not tation of the preS ently available data, may suggest a novel 
. directly observable, but their presence is deduced from their wa y Q f approaching this problem 
5— i ' effect on galactic dynamics, galactic halos are relatively less 
^3 studied structures of galaxies. Various populations, such as 
globular clusters, RR Lyrae, subdwarfs, and other types of 
stars, have been used to trace the halo of the Galaxy, its 
vertic al structure, and its rotation speed (|Kinman et_aL|, In our analy s is we use the seven-year WMAP data 
l 2 004)- Nevertheless, there are still many ambiguities not (iJarosik et all l20ll in the three bands W (94 GHz), V 
only in the main halo constituents, but also in the basic (61 GHz), and Q (41 GHz). Using three WMAP bands 
properties such as, in particular, in rotation. i s important in revealing the possible contribution of the 
The degree to which galactic halos rotate with respect Galactic foregrounds since dust, free- free, and synchrotron 
to the disks is difficult to inve stigate; actually, as sta ted in emission contributes differently in each band. Here we re- 
the most recent study of M31 (jCourteau et all l20Tot). test- mind the reader that the band least contaminated by the 
ing for the rotation of M31's halo is still beyond our reach, synchrotron radiation of the Galaxy is the W-band, which 
also has the highest angular resolution. The CMB map's 



2. The 7-year WMAP analysis 



Send offprint requests to: F. De Paolis, e-mail: general structure in the W-band in the region of M31, 
francesco.depaolis@le.infn.it with the marked 4° radius circle (although our analysis ex- 
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Table 1. Temperature excess in the M31 regions for the non foreground-reduced WMAP maps. 



R, deg, kpc 


Region 


N, pix 


W, T m ± a 


V, T m ± a 


Q, T m ± a 


1.65, 21.4 


Nl + Ml + SI 


324 


63.1 ±5.6 


67.2 ± 5.4 


90.0 ±4.2 




N2 + M2 + S2 


321 


20.3 ±4.7 


17.3 ±4.3 


37.0 ±3.3 




Nl + SI 


205 


45.5 ± 5.7 


38.0 ± 5.3 


64.1 ± 4.0 




N2 + S2 


205 


33.8 ± 5.9 


34.3 ± 5.3 


41.8 ± 4.1 




Ml 


119 


121.4 ± 19.6 


117.6 ± 10.0 


134.3 ± 7.4 




M2 


116 


-7.7 ±7.4 


-12.7 ±6.8 


28.4 ±5.5 


2.40, 31.1 


Nl + Ml + SI 


670 


43.7 ±3.6 


43.5 ±3.4 


66.0 ±2.8 




N2 + M2 + S2 


664 


21.0 ± 3.6 


19.1 ±3.3 


35.6 ±2.7 




Nl + SI 


506 


41.0 ± 3.5 


32.2 ± 3.1 


55.6 ± 2.6 




N2 + S2 


504 


24.5 ± 4.3 


23.6 ± 3.9 


34.8 ± 3.3 




Ml 


164 


73.0 ± 9.4 


78.3 ±9.3 


98.2 ± 7.4 




M2 


160 


11.2 ± 6.3 


5.8 ±6.4 


38.1 ±4.8 


3.20, 41.5 


Nl + Ml + SI 


1176 


36.5 ±2.7 


38.9 ± 2.4 


59.6 ±2.1 




N2 + M2 + S2 


1166 


16.0 ± 2.7 


11.3 ± 2.5 


30.3 ±2.2 




Nl + SI 


980 


35.3 ± 2.7 


34.0 ± 2.3 


53.9 ± 2.1 




N2 + S2 


976 


16.4 ± 3.0 


11.5 ± 2.8 


27.7 ± 2.4 




Ml 


196 


63.8 ± 8.4 


63.2 ± 8.3 


88.1 ± 6.6 




M2 


190 


12.5 ± 5.8 


10.3 ±5.8 


44.0 ±4.4 


4.00, 51.9 


Nl + Ml + SI 


1818 


37.4 ± 2.2 


39.6 ± 2.0 


56.7 ± 1.7 




N2 + M2 + S2 


1808 


1.7 ± 2.3 


-2.4 ±2.1 


16.9 ±1.9 




Nl + SI 


1610 


36.8 ± 2.2 


36.9 ± 1.9 


53.3 ± 1.7 




N2 + S2 


1609 


- 0.4± 2.5 


- 3.8 ± 2.2 


13.5 ± 1.9 




Ml 


208 


64.7.0 ±8.1 


60.8 ± 7.9 


83.6 ±6.4 




M2 


200 


12.5 ± 5.6 


9.4 ±5.6 


43.6 ±4.3 


6.00, 77.8 


Nl + SI 


3748 


29.7 ± 1.5 


27.0 ± 1.4 


44.0 ± 1.2 




N2 + S2 


3749 


11.3 ± 1.7 


7.1 ± 1.5 


25.8 ± 1.3 


8.00, 103.8 


Nl + SI 


6606 


34.3 ± 1.2 


34.7 ± 1.1 


51.5 ± 4.0 




N2 + S2 


6600 


19.2 ± 1.3 


15.0 ± 1.2 


38.7 ± 1.0 



Note s. The radius of the considered annulus is given in degrees and in kpc in the first column; the value of 744 kpc (IVilardell et al.l . 
2010) is adopted for the distance to M31. The considered regions as in Fig. 1 (right panel). The numbers of pixels in each region 
are given. The last three columns show the CMB mean temperatures per pixel of each region in /iK in the W, V, and Q bands, 
respectively, with the corresponding la errors (see text for details). 




Fig. 1. In the left panel, the WMAP W-band towards 
the M31 galaxy. The 8.5° x 8.5° sky field centered at 
(121.17°, -21.57°) with the marked 4° circular region. The 
oblique strip indicates the M31 disk, and the analysis in 
the halo region of M31 galaxy is extended far beyond the 
region indicated in the figure. The detailed geometry (up 
to 8°) used in the analysis is shown in the right panel. 



the major axis is slightly more than about 3° and along 
the minor axis is about 1°. Radio observations have shown 
that the M31 HI disk is more exte nded with resp ect to 
the stellar disk ()Chemin et all . l2009t ICorbelli et all I2010D . 
with a major axis sizes of about 5.6° and a minor axis size 
of about 1.2°. In this paper, the adopted Ml and M2 disk 
regions (Fig. [TJ right panel) have major axis size of 8° and 
minor axis size of 1°; this allows us to retain the warped 
part of the M31 disk in the Ml and M2 regions. Moreover, 
we have checked that it is irrelevant, as far as our analy- 
sis is concerned, to extend the M31 minor axis to 1.2°. The 
mean temperature excess per pixel T mi in /iK/pixel, in each 
region was obtained in each band and is shown in Table 1 
with the la errorQ, along with the number of pixels in each 
area. For convenience, Table 1 gives the temperature excess 
in each M31 region up to 8°, even if our analysis was ex- 
tended to the region around the M31 disk with concentric 
annuli with radii up to 20°. 



tends farther out), is shown in Fig. [T] (left panel). In our 
analysis, we also used the maps provided by the WMAP 
Collaboration with the Galactic disk contribution modeled 
and removed (jGold et al.l . l201ll) . It is always specified in 
the text when we considered these data. To reveal the dif- 
ferent contributions by the M31 disk and halo, the region 
of the sky around the M31 galaxy was divided into several 
concentric circular areas as shown in Fig. Q] (right panel). In 
the optical band the total extent of the M31 galaxy along 



The standard error of the mean is calculated as the standard 
deviation of the excess temperature distribution divided by the 
square root of the pixel number in each region. We have veri- 
fied that, within the errors, the sigma values calculated in that 
way are consistent with those evaluated by using the covariance 
matrix obtained by a best-fitting procedure with a Gaussian to 
the same distribution. 
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3. Results for the M31 disk 

For the M31 disk, our analysis shows that each Ml region 
is always hotter than the corresponding M2 region, as can 
be seen from Table 1. Indeed we find a temperature excess 
contrast (i.e. the difference between the temperature ex- 
cesses per pixel) between the Ml and M2 regions in all three 
WMAP bands that turns out to be about 129±21 /xK/pixel 
within 21.4 kpc (in the W band) and then slightly decreases 
(but remains as large as about 41 ± 10 /iK/pixel at about 
50 kpc). This effect seems to come from the rotation in- 
duced Doppler shift of the gas and dust emission from the 
M31 disk - indeed, the hotter (Ml) region corresponds to 
the side of the M31 disk that rotates towards us. If one 
compares what WMAP data show towards the Ml and M2 
regions with the maps of the M31 thick HI disk ob tained at 
21 cm (jChemin et all . l2009t ICorbelli et all . I2010T ) one sees 
a remarkable superposition of the hot (Ml) and cold (M2) 
regions in both observations. Even if the temperature asym- 
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Fig. 2. The la (green lines), 2a (brown lines), and 3a 
(red lines) excess temperature contrast (in mk/pixel) curves 
(in the W, V, and Q bands) along with the mean profile 
(pink line close to zero) for 500 random control fields (con- 
tinuous lines) and 500 simulated CMB sky maps (dashed 
lines). In red, the observed temperature contrast profile in 
the M31 disk (with la errors) is given. The WMAP maps 
with the Galactic disk contribution modeled and removed 
(foreground-reduced maps) are used here. 



metry between the Ml and the M2 regions looks significant, 
we have to check whether it might be due to a random fluc- 
tuation of the CMB signal. It is indeed well known that the 
CMB sky map has a "patchy" structure characterized by 
the presence of many hot and cold spots with temperature 
excesses up to several tents of /xK/pixel on angular degree 
scales. We therefore considered (Fig. [5] and also the online 
material) 500 control fields and 500 simulated sky maps 
(from the best-fitted cosmological parameters as provided 



2 A detailed study of the frequency dependent temperature 
asymmetry in the CMB arising from different distributions of 
gas and dust in the M31 disk is left to a forthcoming paper. In 
any case, although some inhomogeneity in the disk structure is 
not excluded, there is no reason to assume that it is the sole 



in the WMAP web site) - and evaluated the temperature 
contrast with the same geometry as was used towards M31. 
We also give (red curve) the M31 temperature contrast pro- 
file in the M31 disk. Due to our chosen geometry, each curve 
is given up to 4°. As one can see, the contrast temperature 
profile for the M31 disk is always a nicely smooth curve 
that is close to the 2a curve in the intermediate region of 
about 1.5-2°. Both the control field and the CMB simu- 
lation analyses show that there is a probability of less than 
about 4% that the temperature asymmetry revealed comes 
from a random fluctuation of the CMB signal. Actually, if 
one takes the direction of rotation of the M31 disk into ac- 
count, such a probability reduces (by using the theorem of 
the composite probability) by a factor of two. Finally, we 
mention that we have found that the temperature excess 
contrast of the two M31 disk regions obtained by dividing 
the M1+M2 region with respect to the north-west /south- 
east symmetry axis (the M31 disk axis) turns out to be 
0.008 ±0.012 mK, which seems to further confirm that the 
temperature contrast between the Ml and M2 regions is 
not due to a random fluctuation in the CMB signal. 

4. Results for the M31 halo 

The next step was to enlarge our analysis to the region 
around the M31 disk by considering concentric circular re- 
gions of increasing galactocentric radii (see also Fig. 1, right 
panel). We estimated the difference of the temperature ex- 
cess in the region Nl+Sl in the three WMAP bands with 
respect to that in the region N2+S2. A temperature con- 
trast between the region Nl+Sl with respect to N2+S2 
shows up (see Table 1), and the Nl+Sl region turns out to 
always be hotter than the N2+S2 region. The detected ef- 
fect resembles the one towards the M31 disk, although with 
less temperature asymmetry. In all three bands, the maxi- 
mum temperature contrast reaches a maximum at a galac- 
tocentric distance of about 4° and then decreases slightly. 
It is apparent from the size of the considered regions that 
a contamination of the M31 disk in the regions Nl, N2, SI, 
and S2 can be completely excluded, and also the Galactic 
plane emission cannot account for the observed tempera- 
ture asymmetry since it eventually would make a larger 
contribution towards the upper regions of M31 (while the 
opposite is observed in the data). As for the M31 disk, the 
temperature asymmetry in the M31 halo is indicative of 
a Doppler shift modulated effect possibly induced by the 
rotation of the M31 halo. 

Also in this case we need to check the robustness of our 
results; that is, we have to estimate the probability that the 
temperature asymmetry in the M31 halo is due to a ran- 
dom fluctuation of the CMB signal. In Fig. [3] (see also the 
online material) we have considered 500 control fields and 
500 simulated sky maps (from the best-fitted cosmological 



3 CMB maps were simulated by assuming AT(n) = 
ATcMB(n) ® B{n) + N(h), where ATcmb is a realization of 
the Gaussian CMB field, N(h) is the pixel noise and B(n) is 
the proper beam of the experiment. We made 500 realizations 
of the CMB sky using synf ast routine of HEALPix with the 
best-fit power spectrum as given by the WMAP Collaboration. 
The maps are then convolved with WMAP beams for W, V, 
and Q bands, respectively. Noise realizations (simulated with 
cr = 6.549 mK, a = 3.137 mK, and a = 2.197 mK for W, 
V, and Q-bands, respectively) are added to the beam convolved 
maps in the end. 
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parameters as provided on the WMAP website). As one can 
see, in all three bands, the contrast temperature profile of 
the M31 halo is close to the la curve up to about 10° and 
goes slightly beyond it at about 50 — 60 kpc where the halo 
effect is maximum. This means that there is less than 30% 
probability that the temperature contrast we see towards 
the M3f halo is due to a random fluctuation of the CMB 
signal. We also point out that we have verified that the tern- 





W BANDE 










Fig. 3. The same as Fig. 2 but for the M31 halo (tem- 
perature contrast of the Nl+Sl with respect to the N2+S2 
regions in mK/pixel - red line) for 500 random control fields 
(continuous lines) and 500 simulated sky maps (dashed 
lines). Here the WMAP maps with the Galactic disk con- 
tribution modeled and removed (foreground-reduced maps) 
are used. 



perature asymmetry towards the M31 halo vanishes if the 
adopted geometry is rigidly rotated of an angle larger than 
about I 0° with respect to the M31 major axis, thus giving a 
further indication that the halo temperature contrast effect 
might be genuine and not simply a random fluctuation of 
the CMB. 

We also point out that the use of three WMAP bands 
is useful for revealing the role of the contribution to the 
Galactic foregrounds since each emission mechanism con- 
tributes differently in each band. That the temperature 
contrast seems present in all three bands and is more or 
less the same in each band up to about 10° — 11° indicates 
that the foregrounds are far weaker than the effect. This 
size corresponds to the typical size inferred for the dark 
matter halos around massive galaxies and might open the 
possibility of a new way of studying these systems, galac- 
tic disks, and halos, by using the microwave band. In any 
case, a careful analysis of the Planck data that will be re- 
leased shortly should allow either proving or disproving our 
results. 



5. Discussion and conclusions 

We have shown that a temperature asymmetry in all 
WMAP bands may exist both in the M3I disk and halo in 
the direction of the M3I spin. For the M3I disk, the effect 
is fairly clear, and there is a probability below about 2% 
that it is a random fluctuation of the CMB signal. If real, 



the detected temperature excess asymmetry should be due 
to the foreground emission of the M31 disk modulated by 
the Doppler shift of the disk spin. That the present study is 
really timely is strengthened by considering that the M3I 
galax y is already detected by the Planck observatory (|Add . 
l20Iil )Fl. whereas it did not appear in the WMAP list. These 
are all reasons to expect that the particular effect we discuss 
here can be studied more accurately with Planck data. 

As for the M31 halo, we have shown that, although 
less evident than for the M3I disk, there is some evidence 
of a temperature asymmetry between the Nl+SI and the 
N2+S2 regions that resembles that of a Doppler shift effect 
induced by the M3I halo rotation. We have shown in the 
previous section that there is less than about 30% probabil- 
ity that the detected temperature asymmetry at a galacto- 
centric distance ~ 50 kpc comes from a random fluctuation 
of the CMB signal. 

If one assumes that this temperature asymmetry in the 
M3I halo relies in the M31 itself and is related to the M3I 
halo rotation, one could speculate about the origin of this 
effect. In general, four possibilities may be considered: (a) 
free-free emission, (b) synchrotron emission, (c) Sunyaev- 
Zel'dovich (SZ) effect, and (d) cold gas clouds populating 
the M3I halo 0. To work, the first three effects, assume 
the presence of a rather hot plasma in the halo of M3I. 
Although this hot plasma has not been detected yet, one 
can assume that a certain amount of this plasma can popu- 
late the M3I halo (spiral galaxies are believed to have much 
less hot gas than ellipticals) and may rotate with a certain 
speed. Free- free emission arises from electron-ion scattering 
while synchrotron emission comes mostly from the acceler- 
ation of cosmic-ray electrons in magnetic fields. Both effects 
give rise to a thermal e mission with a rathe r steep depen- 
dence on the frequency (|Bennett et all [2003) that therefore 
should give a rather different temperature contrast in the 
three WMAP bands. The absence of this effect indicates 
that the contribution from possibilities (a) and (b) should 
be negligible. And for (c), even for typical galaxy clusters 
with diffuse gas much hotter than that possibly expected 
in the M3I halo, the rotational scattering effect would pro- 
duce a temperature asymmetry of at most a few /iK/pixel, 
depending on the rotatio nal velocity and the incl ination an- 
gle of the rotation axis dCooravfe Ched . l200l . Actually, 
a possible temperature asymmetry in the CMB data to- 
wards the M31 halo as a consequence of the existence of 
a p opulation of cold gas c louds in its halo was predicted 
in (|De Paolis et all If 99 5a) - possibility (d). Indeed, if the 



4 However, there is no mention of any temperature asymmetry 
in the M31 disk in that paper. 

° We also mention that the numb er and the temperature pro- 
file of radio sources in CMB maps (IGurzadvan et all . I2010f ) ex- 
cludes their significant contribution in the effect under study. 

6 We also considered the possible influence of the ob- 
served high-ve locity clouds, ei t her in the M31 or in our 
galaxy halos (|Westmeier et all 120081 ; iHulsbosch fc Wakkerl 
l!988l : lMorras et all l2000h . bv removing the pixels in the direc- 
tion of each cloud from the analysis. The results obtained do not 
change with respect to those presented here, as expected when 
also considering the relatively low number of pixels involved. 
Also the prop osed e .j ecta by the past interaction of M33 and 
M3I galaxies (Bekki, 2008) cannot play any role in our analysis 
since it would at most have made hotter some pixels in the S2 
region (where the M3I-M33 bridge is located), which is instead 
colder than the SI one. 
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halo of the M31 galaxy contains cold gas clouds, we expect 
them to rotate like the M31 disk (even if, perhaps, more 
slowly), and thus there should be a Doppler shift inducing 
a temperature anisotropy AT between one side of the M31 
halo and the other with respect to the rotation axis per- 
pendicular to the disk. In the case of optically thin halo 
clouds, the Doppler induced temperature anisotropy would 
be AT/T r ~ 2vS f/c , where v is the M31 rotation speed, f 
the averaged cloud optical depth over the frequency range 
{v\ < v < v-i) of a certain detection band, and S the cloud 
filling factor, i.e. the ratio of filled (by clouds) to total pro- 
jected surface in a given field of view. We emphasize that 
the fact that the temperature contrast in Fig. [3] looks ap- 
proximately the same in each band makes a point towards 
either possibility (d) or a random fluctuation of the CMB 
sky (but with a probability, if estimated purely statistically, 
of less of 30% for the last possibility). 

The wealth of data especially in the last decade shows 
that there is good evidence for the presence in the ha- 
los of spiral galaxies of gas in all gaseous ph ases: neu- 
tral, warm atomic, and hot X-ray emitting gas (Brcgman, 
2009). Atomic gas (often identified as HVCs) is observed 
in the radio band (particularly at 21 cm) and through 
absorption lines towards field stars and quasars. The hot 
gas may be detected in X-rays, while searches for cold 
gas clouds in galactic halos are more problematic as are 
searches for t hem by the presence o f a ga mma-ray halo 



(iDixon et all. I 1998t iDe Paolis et all Il999). stellar scin 



tillations ( Moniezl . 120031 : lHabibi et ah L 20101) . obs curation 
events towards the LM C (iDrake fc Cookl. I2003D. o rtho- 
H 2 D+ line at 372 GHz ()Ceccarelli fc Dominikl . I2006D . and 
extreme scattering event s in quasar radio-flux variations 
(| Walker fc War die. 1998) have given no clear indication of 
their presence. 

In conclusion, we showed that our analysis based on 
seven-year WMAP data suggests there is a temperature 
excess asymmetry in the M31 disk is likely due to the M31 
foreground emission modulated by the Doppler shift in- 
duced by the M31 spin. We find that there is less than 
~ 2% probability that the signal up to about 20 kpc conies 
from a random fluctuation in the CMB signal. For the M31 
halo, we also find a temperature excess asymmetry between 
the Nl+Sl and the N2+S2 regions along the expected spin 
direction, suggestive of a rotation induced Doppler shift. 
The effect in the M31 halo is far weaker than for the disk, 
as obviously expected, and more precise data are necessary 
before drawing any firm conclusion. In all cases, this re- 
search may open a new window into the study of galactic 
disks and especially the rotation of galactic halos by using 
the Planck satellite or planned balloon-based experiments. 
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Fig. 4. In the W, V, and Q bands, the la (green lines), 
2a (brown lines), and 3a (red lines) excess temperature 
contrast (in mk/pixcl) curves, along with the mean profile 
(pink line close) for 500 random control fields. The observed 
temperature contrast profile in the M31 disk (with la er- 
rors) is given in red. The non foreground-reduced WMAP 
maps are used here. 




Fig. 5. As above (in the W, V, and Q bands) the la (green 
lines), 2a (brown lines), and 3a (red lines) excess temper- 
ature contrast (in mk/pixel) curves along with the mean 
profile (pink line close) for 500 random control fields. The 
real temperature contrast profile in the M31 halo up to 20° 
(with la errors) is given in red. 



